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ABSTRACT 


Structural  flexibility  effects  exert  a  considerable  influence  upon  the  stabil¬ 
ity  of  booster  orientation  flight  control  systems.  This  paper  presents  a  rela¬ 
tively  simple  method  of  using  root  locus  techniques  to  develop  the  factored 
open  loop  relating  function  for  the  system  including  these  flexibility  effects. 

The  analysis  procedure  is  very  useful  in  the  design  of  these  control  systems 
since  the  solution  to  the  stability  problems  encountered  becomes  evident  from 
the  form  of  the  factored  open  loop  relating  function. 

The  particular  manner  in  which  the  development  takes  place  enables  one 
to  understand  the  effects  of  the  major  design  and  physical  parameters  upon 
the  form  of  the  open  loop  relating  function,  and  in  so  doing,  bridges  the  gap 
between  these  parameters  and  their  effect  on  the  performance  of  the  closed 
loop  system. 

Armed  with  the  basic  understanding  which  the  analysis  procedure  provides, 
it  is  possible  to  interpret  the  performance  of  a  parameter  adjusting  adaptive 
control  system  or  to  formulate  rapidly  the  preliminary  design  of  a  fixed  para¬ 
meter  flight  control  system. 
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SIMPLIFIED  ANALYSIS  OF  FLEXIBLE 
BOOSTER  FLIGHT  CONTROL  SYSTEMS 

by 

Lee  Gregor  Hofmann  and  Allen  Kezer 


INTRCPUCTION 

Structural  flexibility  effects  have  been  extremely  important  in  the  design 
of  flight  control  syatems  for  large  ballistic  missiles  and  boosters.  As  the 
size  of  the  vehicles  has  increased,  the  lightly  damped  structural  modes  of 
vibration  have  moved  into  the  same  frequency  range  that  is  desired  for  the 
flight  control  system  so  that  strong  coupling  effects  have  resulted.  The  nature 
of  these  coupling  effects  has  been  such  that  it  has  been  difficult  to  achieve 
adequate  stability  margins  for  all  flight  conditions.  The  relatively  low  stability 
margins  have  necessitated  complex  and  detailed  analysis  of  the  systems  to 
Insure  that  the  margins  do  indeed  exist  and  that  uncertainties  in  the  knowledge 
of  vehicle  parameters  will  not  lead  to  system  instabilities. 

In  the  process  of  the  development  of  an  adaptive  control  system  to  over¬ 
come  some  of  these  obstacles,  it  was  necessary  to  develop  simplified  analysis 
techniques  which  would  lead  to  a  better  understanding  of  the  effects  of  parameters 
upon  the  performance  of  the  system.  These  simplified  analysis  techniques  are 
usefxil  not  only  in  the  design  of  adaptive  control  systems  but  also  in  the  prelim¬ 
inary  design  of  a  basic  control  system  to  which  the  adaptive  features  may  be 
added  if  necessary. 

One  of  the  prime  considerations  in  the  design  of  the  flight  control  systems 
is  the  fact  that  the  sensors  which  provide  the  means  of  closing  the  feedback 
loops  (gyros,  accelerometers,  etc  J  sense  the  local  bending  at  the  sensor 
station  as  well  as  the  motion  of  the  ideal  rigid  body.  In  some  cases,  a  very 
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significant  improvement  can  be  made  in  the  closed  loop  bending  dynamics  by 
use  of  proper  feedback  of  the  bending  information  that  is  sensed  by  these  instru¬ 
ments.  This  can  result  in  a  very  significant  reduction  in  the  bending  motion 
and  hence  in  the  dynamic  loads  which  are  applied  to  the  structure. 

In  any  feedback  control  system,  the  feedback  signal  must  have  the  proper 
phase  and  amplitude.  Usual  root  locus  or  frequency  response  methods  of  analy¬ 
sis  which  are  used  to  insure  the  existence  of  the  proper  phase  and  amplitude 
require  that  the  open  loop  poles  and  zeros  of  the  system  transfer  function  be 
known  in  factored  form.  This  paper  employs  a  relatively  simple  and  systematic 
method  of  obtaining  the  factored  open  loop  transfer  function  for  a  booster  con¬ 
trol  system  using  simplified  equations  of  motion  and  a  root  locus  factorization 
technique. 

The  development  illustrates  the  application  of  this  technique  to  the  case 
of  a  system  which  controls  the  orientation  of  the  undeflected  centerline  of  the 
booster.  In  this  system  the  feedback  quantities  are  angular  displacement  and 
rate  of  angular  displacement  with  respect  to  inertial  space  of  the  flexed  center- 
line  of  the  booster  at  one  or  more  locations  along  the  centerline. 

The  feedback  of  other  quantities,  such  as  normal  acceleration,  may  be 
included  in  the  control  system  by  developing  an  applicable  set  of  equations 
using  the  method  of  derivation  employed  in  this  paper.  The  simplifying 
assumptions  which  make  this  systematic  analysis  possible  do  not  cause  any 
significant  loss  in  accuracy  in  the  description  of  the  physical  booster,  and 
because  of  the  dynamic  similarity  of  the  booster  modes  of  motion,  the  analy¬ 
sis  is  effective  in  reducing  the  apparent  complexity  of  the  problem.  The 
effects  of  the  important  booster  characteristics  and  design  parameters  upon 
system  performance  are  quite  evident,  and  parametric  studies  are  greatly 
facilitated  by  this  particular  method  of  obtaining  the  open  loop  relating  function. 

ANALYSIS  OF  THE  OPEN  LOOP  PERFORMANCE  FUNCTION 
OF  THE  FLIGHT  CONTROL  SYSTEM 

In  order  to  analyze  the  effects  of  parameter  variations  on  the  closed  loop 
performance  of  the  booster  flight  control  system,  the  open  loop  relating 
function  must  be  known  in  factored  form.  To  proceed  toward  this  end,  a 
mathematical  description  of  the  control  system  must  be  formulated  from  a 
functional  description  of  the  system. 
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Functional  Description  of  Control  System 

A  functional  block  diagram  of  an  orientation  flight  control  system  is  presented 
in  Fig.  1.  This  system  utilizes  feedback  of  angular  orientation  information  as 
sensed  by  rate  and  rate  integrating  gyros.  In  many  cases,  the  orientation 
measurement  may  be  performed  by  the  inertial  measurement  unit  of  the  guidance 
system  instead  of  a  separate  rate  integrating  gyro.  Angular  orientation  or 
angular  rate  may  be  commanded  to  the  system  by  suitable  usage  of  the  orien¬ 
tation  measuring  instrument.  The  input  and  output  variables  of  the  system 
are  different  for  the  different  commands,  but  the  loop  dynamics  and  analysis 
procedures  are  unaffected. 

The  rate  and  rate  integrating  gyros  may  be  located  at  different  stations 
along  the  booster  body.  If  necessary,  the  outputs  of  several  rate  (or  rate  inte¬ 
grating)  gyros  may  be  summed  in  such  a  way  that  the  rigid  body  component  is 
the  same  as  that  sensed  by  a  single  gyro,  while  the  effective  bending  mode 
slope  "sensed"  for  each  mode  may  be  changed  by  suitable  adjustment  of  the 
contribution  made  by  each  gyro. 

The  introduction  of  compensation  into  the  control  loops  provides  several 
additional  system  variables  which  in  turn  make  it  possible  to  use  different 
approaches  for  stabilizing  the  different  modes  considered  in  the  design.  For 
generality,  compensation  units  are  placed  at  three  locations  in  the  two  control 
loops  of  Fig.  1. 

In  this  study,  it  is  assumed  that  a  gimballed  nozzle  is  used  for  the  thrust 
vector  control  and  that  a  hydraulic  actuator  is  used  to  provide  the  required 
nozzle  deflection.  For  clarity,  a  first  order  lag  approximation  is  used  for 
the  hydraulic  actuator  and  a  second  order  system  for  the  nozzle  dynamics. 
Accelerations  of  the  nozzle  gimbal  point  which  produce  reaction  torques  on  the 
nozzle  and  couple  the  rigid  body  and  bending  modes  with  the  nozzle  motion  are 
neglected,  but  may  be  included  by  mathematical  operations  which  are  similar 
to  those  developed  in  this  paper. 

Mathematical  Description  of  System  Dynamics 

The  equations  of  motion  which  have  been  developed  for  the  booster  air¬ 
frame  are  quite  complex*  and  use  of  these  complete  equations  is  cumbersome. 

To  facilitate  analysis  and  interpretation  of  parametric  effects,  the  equations 
may  be  simplified  without  any  significant  loss  in  accuracy.  In  Appendix  A, 

*  See  Appendix  A 
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INTEGRATING 

GYRO 

COMPENSATION 

(>g) 

(C3) 

COMPENSATION 

(c,) 


ACTUATOR 

(o) 


BOOSTER 


(■) 


1  MOTION  IITH  RESPECT  TO  INERTIAL  SPACE  SENSED  BY  THE  INTEGRATING  GYRO  (^,,) 

2  MOTION  OF  BOOSTER  UNDEFLECTED  CENTER  UNE  HTH  RESPECT  TO  INERTIAL  SPACE  (0) 

3  MOTION  NITM  RESPECT  TO  INERTIAL  SPACE  SENSED  BY  THE  RATE  GYRO  [4,g) 

4  INERTIAL  LOADING  EFFECT  UPON  ACTUATOR 


Fig.  1.  Functional  block  diagram  of  boosftr  flight  control  systtm. 


NOZZLE  DYNMMCS 


THE  EQUATIONS  FOR  THE  ABOVE  RELATING  FUNCTIONS  ARE  LISTED  BELOW 


Fig.  2.  Block  diagram  roprosontation  of  tho  timplifiod  boostor  dynamics. 
(Tho  offocts  of  throo  bonding  modot  oro  includod.) 
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Fig.  3.  MotliMiatical  block  diagram  of  tho  opon  loop  booslof  control  systom  for  on  arbitrary 
number  of  bonding  modps  -  simplifiod  oquotions  of  motion. 


the  complete  linearized  equations  have  been  simplified  to  a  degree  which  permits 
the  relating  function  for  the  output  of  the  Booster  Dynamics  functional  block 
of  Fig.  1  to  be  represented  by  the  block  diagram  of  Fig.  2  when  the  (rg)  notation 
is  substituted  in  the  subparenthesis  of  9^  ^  and  y  The  relating  function 


notation.  This  diagram  in  Fig.  2  is  a  general  representation  of  the  booster 
equations  from  a  nozzle  deflection  (6)  to  the  motion  sensed  by  a  single  gyro 


The  relating  function  paths  of  the  rate  damping  and  orientation  control 
loops  can  be  combined  into  one  path  for  each  booster  mode.  The  mathematical 
block  diagram  of  the  resulting  open  loop  relating  function  is  shown  in  Fig.  3. 

The  opening  in  the  control  loops  depicted  in  Fig.  3  corresponds  to  an  opening 
ahead  of  the  sensitivity  adjustment.  in  the  forward  portion  of  the  control 

loops  of  Fig.  1. 

The  parallel  paths  shown  in  Figs.  2  and  3  result  from  the  fact  that  the  gyros 
sense  a  component  of  angular  motion  due  to  the  local  slope  of  each  bending  mode 
as  well  as  the  angular  motion  of  the  rigid  body.  The  factor  which  represents 
the  local  slope  of  each  bending  mode  (X^  ^  )  appears  as  a  separate  sensitivity 
factor  in  the  mathematical  block  diagrams^of  Figs.  2  and  3.  The  slope  of  the 
i^^  bending  mode  at  the  integrating  gyro  station  is  denoted  by  Xj^g,  and  at  the 
rate  gyro  station  by  X  In  Fig.  3  the  parallel  paths  also  contain  the  compen- 
sation  elements  which  are  unique  to  either  the  rate  damping  loop  or  the  orientation 
control  loop. 

The  open  loop  performance  function  of  the  booster  control  system  may  be 
written  by  inspection  of  Fig.  3: 

(PF)b«.  =(RF)b[8;8;]  Soc(PF)c2  [8';8,]  (Pn„[8,;8] 

epan 

loop 

The  parallel  paths  of  the  open  loop  relating  function  which  compose 
(RF)  giij  and  represent  the  components  of  the  control  system  feedback 
signal  due  to  nozzle  deflection,  must  be  reduced  to  a  single  path.  Algebraic 
reduction  of  the  expressions  for  these  parallel  paths  to  an  expression  for  a 
single  path  results  in  a  ratio  of  polynomials  with  a  factored  denominator  and 
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an  unfactored  numerator.  The  numerator  polynomial  must  be  factored  before 
proceeding  with  an  analysis  of  the  closed  loop  control  system. 

The  development  of  the  expressions  for  obtaining  the  factors  of  the 
numerator  polynomial  is  based  on  the  algebraic  manipulation  of  the  sum  of 
two  polynomials  into  a  familiar  unity  feedback  expression  to  which  root  locus 
factoring  techniques  may  be  applied. 

For  example: 


A,  B,  C,  D,  q^,  and  q^^^  are  functions  of  p;  A,  B,  C,  and  D  are  factored 
polynomials  in  p. 

huL  -  A  xC  .ADJ_CB 
"B^D*  BD 


If  A,  B,  C,  and  D  are  factored  polynomials,  the  quantity  of  the  inner  brackets 
(  ]  may  be  factored  by  a  root  locus  where  the  "open  loop"  poles  are  the  factors 
of  A  and  D  and  the  "open  loop"  zeros  are  the  factors  of  B  and  C.  Because  of  the 
minus  unity  exponent  on  the  term  in  the  outer  bracket  |  | ;  the  locus  of  the 
poles  of  the  root  locus  factorization  becomes  a  locus  of  the  numerator  zeros 
of  the  qout/Qjij  transfer  fimction. 
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The  determination  of  the  zeros  of  the  relating  function  of  the  open  loop 
booster  control  system  shown  in  Fig.  3  involves  the  successive  application 
of  this  technique  as  the  parallel  path  for  each  bending  mode  is  added  to  the 
analysis.  The  equations  for  the  determination  of  the  zeros  are  presented  in 
Equation  Summary  1  for  the  case  in  which  the  compensation  imits  c^  and  c^ 
have  transfer  functions  of  unity.  *  The  expressions  for  including  one  bending 
mode  in  parsillel  with  the  rigid  body  mode  are  given  by  equations  l-3(a).  (b). 
(c),  and  1-4. 

The  extension  of  this  result  to  allow  the  addition  of  the  effects  of  a 
bending  mode  when  an  arbitrary  number  of  bending  modes  have  already  been 
included  in  parallel  with  the  rigid  body  mode  is  given  by  Eqs.  l-7(a),  (b),  (c), 
and  1-8. 


*The  results  of  a  similar  development  for  a  more  general  case  in  which  the 
compensation  units  c^  and  C2  do  not  have  performance  functions  of  unity  is  pre¬ 
sented  in  Appendix  B.  However,  the  results  there  are  mainly  of  academic  interest, 
since  the  effects  of  these  compensation  unit  dynamics  may  be  accurately  repre¬ 
sented  by  simpler  means.  The  simpler  method  involves  defining  and  calculating 
effective  values  of  the  local  bending  mode  slopes  "sensed"  by  the  rate  and  inte¬ 
grating  gyros.  These  effective  values  of  the  slopes  are  equivalent,  in  the  system 
performance  sense,  to  the  combination  of  the  actual  values  of  the  locsd  bending 
mode  slopes  sensed  by  the  gyros  and  the  effects  of  the  compensation. 
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Asauming  that  the  perfermance  functions  of  the  cj  and  C3  compensation  units  ate  equal  to  unity,  the  fol¬ 
lowing  espreasion  for  (RF)^  [  $;$*.]  can  be  developed  ftom  the  mathematical  block  diagram  of  Fig.  3  and  the 
relating  functions  listed  in  Fig.  2 


\  *1"1  0,  ) 

The  first  term  in  the  above  eipression  is  the  contribution  of  the  rigid  body  motion  to  the  ig.  signal;  the 
second  term  is  the  contribution  of  the  lat  mode  motion.  Placing  the  above  eipression  over  a  common  denom¬ 


inator  gives 


V  "1  TS.  J  /  ^ 


Introducing  s  simplifying  notation 


Ht5;5] 


[A,(p)+K,  B,{p)] 


Equation  Summary  1  (Page  1  of  4) 

[/cvclopment  of  the  Flexible  Booster  Relating  Function  — 
Without  Compensation  ((PF)jj  ={PF)jj  =1] 


_  ~^ig]  ^b[S;q^] 
’  Sk[S;§] 


(l-3b) 


(1-3C) 


In  Eqa.(l*2)  and  (1-3)  che  five  zeroa  of  (he  relatiajt  function  are  contained  in  the  bracketed  term.  One  zero 
reaulta  from  (he  rate  gyro  “differentation**;  two  from  placing  the  polynomials  over  a  common  denominator,  and 
two  ''(ail-waga*dog*’  zeroa  from  (he  inertial  reactions  of  the  gimballed  nozzle.  Care  must  be  taken  not  to 
confuse  (he  "tail-wags-dog”  zero  of  the  booster  which  is  a  factor  of  (R^)b[S;S*]  *‘*b  the  "tail-wags-dog” 
zero  of  the  ith  mode  which  is  a  factor  of  (RF)^[  • 

By  rearranging  the  bracketed  term  of  Eq.  (1-3).  the  polynomial  summation  A|  +1^]  Bj,  may  be  put  into  the 
familiar  unity  feedback  form,  multiplied  by  a  factored  polynomial,  l/K^ 


(Rf')b[8;s;] 


Sb[8;§] 

®i  1 

''*7 

1 

K,B, 

•1 


(1-4) 


The  quantity  in  (he  inner  brackets  of  Eq.  (1*4)  is  the  unity  feedback  expression  which  is  readily  factored 
by  root  locus  techniques. 

The  final  factored  form  of  the  relating  function  is  obtained  by  completing  the  other  algebraic  operations 
that  are  indicated  in  Eq.  (1-4);  noting  especially  the  minus  unity  exponent  of  the  outer  brscketed  term. 


The  effects  of  additional  bending  modes  can  be  included  in  the  relating  function  by  similar  operations. 


In  adding  the  second  mode  the  equations  become 


for  first 
and  saeend 
modes 


=(RF)b[8;8*) 


for  first 
mode 


Equation  Summary  1  (Page  2  of  4) 

Development  of  the  Flexible  Booster  Relatiiy  Function  — 
Without  Compensation  “  (PF)^^  ®1]  ■ 


o 
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Ezainiauion  of  the  above  equation  suggeata  that  it  ia  poaaible  to  generalize  it  immediately  to  the  caae  where 
the  nth  mode  ia  being  added  to  the  ezpreaaion  for  the  (n>  1)  mode  relating  function.  Additionally,  it  can  be 
Been  that  it  ia  not.  neceaaarily  the  mode  n  which  muat  be  added  at  thia  point  aince  any  of  the  other  (n*  1)  modea 
might  have  been  uaed  in  ita  place.  Therefore,  the  mode  added  will  be  denoted  j  to  indicate  that  the  modea  need 
not  be  included  in  order  according  to  their  mode  index. 

The  method  for  the  calculation  of  the  relating  function  which  includea  n  bending  modea  when  the  relating 
function  which  includea  (n  •  1)  modea  ia  known,  is  derived  in  the  following  operationa. 


{RF)b[ 


^l«j  ^b[S;qj]  0  2  )  » 

x-s'l  =(RF)kf  - yr - 5 — ^  (  '  * P  j  (l-5b) 

»;»eJ  forn  J  fe,  („.i)  /  ,  ^  ^  _pL  \  '  ^Igi  ' 

mod..  mod..  2)  ' 

I  )  0:  <U; 


The  firat  term  in  the  above  expreaaion  is  the  contribution  of  the  rigid  body  motion  and  (n*  1)  of  the  n  bending 
modes  to  the  S'  signal;  the  second  term  is  the  contribution  of  the  jth  mode  motion. 

Substituting  the  factored  form  of  (RF)b[$.5*]  modo.*"*^  P**C'“*  ‘I**  ‘hose  expression  over  a 

common  denominator  gives 

.mm.  - 

(RF) b C  S;«; ]  , . 

I  •  >  lit  ^ 


f/ the  (2n  +  1)  sore,  of  \  /  2Cj  p*  \ 

(RF)k  5.5.]  h,  (n  •  1)  mod..  )  \  "i  **  0!  J  f 


•1  ' 


'  V  "I'l 


Equation  Summary  1  (Page  3  of  4) 
Development  of  the  Flexible  Booster  Relsting  Function  — 
Without  Compensation  I(PF)g^ 


latroduciai  •inplifyini  aocuioa 


(RP)k  [*;«']  for 


Ttr 


— r-[Aj  (p)+Kj  Bj  (p)] 

}  "  /,  P*  \  ’ 

•4ai  p*  tl  (I  —  p  +  ^  ^  1 


'  II  ®i  "j 

X  [rti#  (Jn  ♦  1)  i«ro»»f(RF)^[5.jj  ]  for  (n-l)  nedat] 


Bj(p)  =p2  (l  )(  1  *(SR)^  j^p  )  y  (l  p  ) 

«K,.  >9f  '  <  '  '  '  0|  a)|  ' 


(1-7) 


(l-7«) 


(l-7b) 


*  ■  hare] 


The  2n  +  1  ceroe  of  Eqe.  (1*6)  and  (l-7«)  mnet  be  eipieeeed  i»  the  foim 


H[S;g] 


forn  •  "  /.  2{i  d2  \ 


' 

'^iTT 

1 

B: 

'  ''i  ®i 

'•"iTsf 

. 

or(1  *rp). 


-1 


(l-7c) 


(1-8) 


The  Older  of  the  poljmomiaU  Aj  (p)  and  Bj(p)  is  (2i»  +  3)  for  n  =2,  3,  4,  ... .  The  eiplasacioii  of  the  oripia  of 
the  aeros  is  the  same  as  before,  bnc  in  addition,  two  aeios  appear  for  each  bending  node  added. 


Eqnation  Summary  1  (Page  4  of  4) 
Development  of  the  FlesiUe  Booster  Relating  Faaction  — 
Without  Compensation  l{PF)j^  ”(**^)e3 


13 


Example  Illustrating  the  Method  of  Applying  the  Simplified  Analygis 

The  usefulness  of  this  simplified  analysis  in  qualitative  applications  cannot  be 
fully  appreciated  until  one  witnesses  a  demonstration  of  its  power  and  simplicity 
in  an  actual  example.  It  is  the  experience  of  the  authors  and  our  associates  that 
the  routine  of  formal  operations  indicated  in  equations  1-4  and  1-8  of  Equation 
Summary  1  are  easily  carried  out  to  obtain  the  factors  of  the  numerator  polynomial 
of 


Problem:  Find  the  factors  of  the  numerator  polynomial  of  ® 

case  in  which  the  first  and  second  bending  mode  effects  are  included 
in  the  analysis. 


Situation:  1.  The  orientation  and  rate  damping  control  loops  contain  the  same 

compensation  which  is  lumped  into  the  C2  element  of  Fig.  1.  The  c^ 
and  Cj  compensation  transfer  functions  are  therefore  equal  to  unity. 


2.  The  damping  of  the  bending  modes  is  assumed  to  be  negligible. 

£,=52^0 

3.  The  "tail-wags-dk>g"  frequencies,  of  the  rigid  body, 

0  12 

first  bending,  and  second  bending  mode,  respectively,  are  related 
by  the  following  inequality.  * 

4.  The  first  and  second  bending  mode  frequencies  are  such  that 


*  It  can  be  shown  that  the  "tail-wags-dog"  frequencies  are  almost  invariably 
related  to  one  another  by  ihe  above  inequality.  The  greatest  of  these  frequencies 

is  slightly  less  thaiW  —  - .  Differences  between  the  frequencies  in  the 

inequality  are  small. 

An  interesting  fact  is  that  the  "tail-wags-dog"  frequencies  are  relatively 

constant  during  a  constant  thrust  flight  profile  as  opposed  to  the  bending  mode 

frequencies  which  may  increase  over  100%  during  the  flight  profile. 

*  th 

a^  is  the  undamped  natural  frequency  of  the  i  bending  mode 

when  the  rocket  engine  is  thrusting;  this  is  also  the  resonant  frequency  when  the 
damping  of  the  i^^  mode  is  negligible. 
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^2  ^  ^*0 
The  first  bending  mode  resonant  frequency  is 
less  than  the  lowest  "tail-wags-dog"  frequency, 
and  the  second  bending  mode  frequency  is  greater 
than  the  highest  "tail-wags-dog"  frequency. 

5.  The  sign  of  the  ratio  of  the  bending  mode  slope 
sensed  by  the  gyros  at  their  respective  stations 
to  the  bending  mode  displacement  at  the  nozzle 
gimbal  station  is  given  for  each  mode  by: 

o)  ^•92^'^«2 

b)  '*  PowHve 

c)  “^rg/*^gi 
Qualitative  Solution: 

The  zeros  of  including  the  first  bending  mode 

in  parallel  with  the  rigid  body  mode  are  determined  by  one 
root  locus  operation.  The  effects  of  the  second  bending 
mode  in  parallel  with  the  combination  of  the  rigid  body  and 
first  bending  mode  are  added  by  a  second  root  locus  oper¬ 
ation.  The  steps  of  the  qualitative  solution  of  the  problem 
are  outlined  below.  These  correspond  to  the  mathematical 
operations  indicated  in  Eqs.  1-4  and  1-8  of  Equation  Sum¬ 
mary  1 .  Figures  4  (a-d)  are  graphical  representations  of 
these  steps. 

1.  Since  (PF)  =  (PF)  =  1,  the  equations  of  Equation 

Summary  1  are  applicable  to  this  problem.  The  first 
bending  mode  will  be  added  first.  Plot  the  open  loop 
roots  of  the  term  in  the  inner  brackets  of  Eq.  1-4, 

B^(p)  /  Aj^(p)  on  the  complex  plane.  The  poles  are 
the  roots  from  Eq.  l-3(a)  and  the  zeros  are  the  roots 

fromEq.  l-3(b).  Since  /d  =X  /A  , 

*1  *1  *^“1  »1 

A^(p)  and  B^ip)  contain  common  factors  of  [  l+piSR)^^] 
which  divide  out. 


Fig.  4(a).  Root  locus  dolarminolion  of  tho  floxibis  boostor  relating  (unction  -  steps  1  through  3. 


RESULT  IS  (RF)b{«;a''1  INCLUDING  THE  FIRST  BENDING  MODE  ONLY 


Fig.  4(b).  Root  locus  determination  of  the  flexible  booster  relating  (unction  -  steps  4  and  5. 
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2.  Find  the  roots  of|Kj(Bj  /  A^)  /  (  1+Kj  (B^  /  A^)  ]| 
by  the  root  locus  technique.  Care  must  be  taken 
to  choose  the  proper  angle  criterion  for  the 
locus.  The  locus  criterion  for  this  part  of  the 
problem  is  180°.  The  root  locus  is  easily 
sketched  by  inspection  since  all  of  the  open 
loop  roots  lie  on  the  imagin-Ary  axis  in  the  "p" 
plane.  Note  that  increasing  the  magnitude  of 

=  A.  corresponds  to  increasing  the  magni- 


rg 


1 


ig 


1 


tude  of  the  open  loop  gain  and  therefore  causes 
the  closed  loop  poles  to  move  away  from  the 
open  loop  poles  along  the  imaginary  axis  toward 
the  zeros. 


3.  After  finding  the  roots  of|Kj(Bj/A^)  /  (  l+K^CB^/A^)  ]|, 
multiply  through  by  the  remaining  term  in  the 
outer  brackets,  l/K^B^.  Cancellations  with  all 
open  loop  zeros  occur,  and  the  sensitivity  of  the 
outer  bracketed  term  becomes  unity. 


4.  Invert  the  bracketed  function.  The  minus  unity 

exponent  on  the  brackets  corresponds  to  changing 
the  poles  remaining  after  step  3  to  zeros.  These 
zeros  are  the  zeros  or  factors  of  the  numerator 
of  the  relating  function  (RF)j^j  g.  including  the 
first  bending  mode.  ^ 

5.  Multiply  the  function  resulting  from  step  4  by  the 
sensitivity  and  characteristic  poles  outside  of  the 
outer  bracketed  function  (add  these  characteristic 
poles  to  the  root  locus  plot).  The  resulting  plot 

is  now  (RF'jjj  g,  giij  including  the  first  bending  mode. 

To  add  the  effects  of  the  second  mode,  apply  the  operations  of  Eq.  1-8 
for  n=j  =2. 
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RESULT  IS  (RF)h[a.a^»]  INCLUDING  THE  FIRST  AND  SECOND  BENDING  MODES. 


Fig.  4(d).  Root  locus  dolorminotion  of  tho  floxiblo  boostor  roloting  function  -  stops  9  and  10. 
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6.  Proceed  as  in  step  1  this  time  plotting  Bg/Ag 
as  the  open  loop  roots  of  the  term  in  the  inner 
brackets.  Note  in  Eq.  l-7(a)  that  the  roots 
of  Ag  are  the  2n+l(n*2)  poles  remaining  after 
the  operations  of  step  3  plus  the  two  roots 
which  are  the  characteristic  of  the  bending 
mode  being  added  (second  bending  mode). 

In  Eq.  l-7(b),  B2  has  roots  as  follows:  2  at 
the  origin,  2  at  the  "tail-wags-dog"  fre¬ 
quency  of  the  bending  mode  being  added, 

2n-2(n=2)  corresponding  to  the  characteristic 
roots  of  the  bending  modes  already  included, 
and  one  from  the  effect  of  different  rate  gyro 
and  rate  integrating  gyro  pickups  of  the 
second  bending  mode  slope  component  of 

angular  deflection.  Since  X.  X  , 

ig2  rg2 

(  1  +p(SR)^^  Vg2^^^S2^  '  1  +p(SR)^dl 

of  Ag  do  not  divide  out  as  did  the  corresponding 
factors  for  the  first  bending  mode.  Steps  7 
through  10  are  similar  to  steps  2  through  5. 

In  step  10,  however,  an  additional  pair  of  char¬ 
acteristic  poles  appear  outside  of  the  outer  i 
bracketed  term.  These  are  the  characteristic 
poles  of  the  second  bending  mode.  The 
resulting  plot  is  now  6-6") 

the  effects  of  the  first  and  second  bending  modes. 

Following  are  several  important  observations  to  note  as  a  consequence 
of  this  example;  some  are  directly  obvious,  the  others  result  from  a 
bit  of  qualitative  experimentation  with  the  parameters  of  the  system. 

1 .  The  open  loop  singularities  used  in  a  locus  operation  are  either 
the  singularities  of  the  relating  functions  for  the  booster  modes 
listed  in  Fig.  2.  or  the  roots  of  the  numerator  polynomial  of 
(RF)^j  g.  g  iij  which  were  found  as  a  result  of  the  previous  locus 

operation. 
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2.  The  form  of  the  root  locus  operation  is  such  that  the  zeros  due 

to  the  inclusion  of  the  i^^  bending  mode  proceed  from  the  open 

th 

loop  pole  locations  of  the  i  bending  mode  with  increasing  mag¬ 
nitude  of  the  i  mode  slope  sensed  by  the  rate  gyro,  A  .  The 

rgi 

parameter  is  a  design  variable  because  its  value  is  deter¬ 
mined  by  the  placement  of  the  rate  gyro  along  the  booster  body. 

th 

3.  The  direction  in  which  the  zeros  due  to  the  inclusion  of  the  i  " 

mode  proceed  from  the  i*^  mode  pole  locations  is  determined 

f  h 

by  the  sign  of  .  the  relative  value  of  the  i  bending  mode 
resonant  frequencV  to  the  i^*^  mode  "tail -wags -dog"  frequency,* 
and  the  value  of  A.  /  A  .  **  The  parameter  A.  is  a  design 

igi  rg^  igj^ 

variable  because  its  value  is  determined  by  the  placement  of 
the  rate  integrating  gyro  along  the  booster  body. 


4.  If  A 


rg.  \g.  ^G. 


for  all  the  modes  included  in  the  analysis,  one 


zero  of  occurs  on  the  real  axis  at  -l/(SR)  The 

rest  of  the  zeros  occur  in  conjugate  imaginary  pairs,  in  real 
pairs  symmetric  about  the  origin,  or  in  complex  conjugate 
foursomes  which  are  symmetric  about  the  imaginary  axis. 


5.  If  A 


=  A.  =  for  all  the  modes  included  in  the  analysis,  the 
rgi  igi  G. 

phase  of  /  1 1  at  each  resonant  frequency 

can  be  easily  calculated  knowing  only  the  parameters  of  the  modes 
and  the  mode  slope  sensed  by  the  gyros,  Aq  .  In  a  similar 
fashion  the  phase  or  range  of  phase  at  each  (>f  the  zeros  is  easily 
calculated  even  though  the  positions  of  the  zeros  are  only  quanta- 
tively  known. 


*i.  e.  w.Wa.  less  than  or  greater  than  u  . 

1  v  1 

**Other  parameters  assumed  constant. 
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The  assumption  that  X 


rg 


i  '  % 


=  X„  is  often  useful  to  make  even  if 
Gi 


the  statement  is  not  strictly  true  because  of  the  substantial  "computational" 
simplification  which  is  introduced.  This  assumption  does  not  effect  the 
position  of  the  zeros  relative  to  the  poles  in  the  direction  of  the  imaginary 
coordinate.  The  parameter  ratio  X.  /  ^  p  controls  the  lateral  location, 

or  real  coordinate  of  the  zeros  with  respect  to  the  poles.  This  control  is 
usually  rather  weak  except  when  1  /  (SR)  tends  to  zero  and  /or  jx^^^  /  X^^ 

tends  to  be  large  with  respect  to  unity.  When  the  control  of  X,_  /  X  ^  is 

igi  rg. 

weak,  its  effect  does  not  alter  the  basic  nature  of  the  problem,  and  therefore 
the  assumption  is  useful  when  only  qualitative  understanding  of  a  problem  is 
required. 

The  statement  X  =  X. _  is  strictly  true  when  both  rate  and  rate  inte- 
rgi  igj 

grating  gyros  are  placed  at  the  same  station  along  the  booster  body. 

A  figure  which  emphasizes  points  4  and  5  above,  follows.  This  will  show 
the  zero-pole  configurations  that  6fi")  assume  for  all  possible 

combinations  of  X^  /0^  greater  or  less  than  zero  and  u^iJS^/u^  greater  or 

less  than  unity’*'  when  one  bending  mode  is  included  in  the  analysis.  Four 


.th 


’‘‘If  the  mode  slope  parameters  X  ^  and  X.  or  are  divided  by  the  i 

rgj  ig^  G. 

mode  displacement  at  the  gimbal  station  if)  ,  the  resulting  ratios  are  independent 

*i 

of  the  manner  in  which  the  mode  shapes  are  normalized.  It  can  be  shown  that 
the  sign  of  this  ratio  indicates  the  phase  of  the  r"  mode  component  of  motmn 
sensed  by  the  gyro  relative  to  the  phase  of  the  force  being  applied  to  the  i'^  mode. 

The  expression  greater  or  less  than  w  can  be  non-dimensionalized 

1  ™  1 

by  dividing  through  by  u  to  produce  a  neater  mathematical  statement.  This 

^i 

ratio  indicates  the  dominance  of  the  transverse  component  of  thrust  acting  on 
the  ith  mode  over  the  transverse  component  of  the  inertial  reaction  force  due 
to  gimballing  the  nozzle  at  the  resonant  frequency  of  the  i*"  mode  when 
^“i'V^  /  )<1.  When(wj^a^  /  )>  1  this  ratio  indicates  that  the  transverse 

inertial  reaction  force  is  dominant  over  the  transv^erse  component  of  thrust  acting 
on  the  i*"  mode  at  the  resonant  frequency  of  the  mode.  Note  that  the  transverse 
inertial  reaction  force  on  the  i^“  mode  is  180*  out  of  phase  with  the  transverse 
component  of  the  thrust  acting  on  the  iih  mode  when  the  nozzle  is  gimballed  at  an 
arbitrary  frequency. 
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basic  combinations  are  possible;  these  are  shown  in  Fig.  5(a*d).  Alternate 
features  which  can  arise  in  cases  3  and  4  of  Figs.  5(a-d)  are  shown  if  Figs. 
5(e-g).  The  method  by  which  these  figures  were  constructed  exactly  parallels 
steps  1  through  5  of  the  previous  example  which  showed  the  method  for  obtaining 
the  zeros  of 

c 

Table  1  shows  the  phase  or  range  of  phase  at  each  singularity  of 
(RF)j^l  g .  5"j  /  [  l+p(SR)  for  the  different  cases  based  on  observations  of 

Fig.  5(a-g). 

The  areas  of  the  complex  plane  in  which  the  zeros  are  found  for  the  dif* 
ferent  cases  are  compiled  in  Table  2.  This  table  is  also  based  on  observations 
of  Fig.  5(a-g). 


The  significance  of  tables  such  as  1  and  2  greatly  increases  when  they  are 
constructed  for  cases  which  include  more  than  one  bending  mode  in  the  analysis. 
When  several  modes  are  considered,  it  will  be  apparent  to  the  control  system 
designer  that  this  classification  procedure  is  helpful  in  keeping  track  of  such 
things  as  the  possible  pole-zero  configurations  and  phase-at-the-singularities 


of  (RF)  |j( 

When  all  possible  combinations  of  the  relative  values  of  Xq  to  zero 
and  of  to  unity  for  all  the  modes  included  in  an  analysis  are  quali¬ 

tatively  investigated,  patterns  of  zero  locations  relative  to  the  poles  are  evident. 
These  patterns,  which  are  uniquely  specified  in  terms  of  the  relative  values  of 
these  parameters,  can  be  translated  into  tables  similar  to  1  and  2.  Subsequent 
use  of  the  tables  requires  knowledge  of  Xq  /  0^  and  ^  “  z  modes 

of  the  particular  booster  of  interest,  but  does  not  require  even  a  qualitative 
factorization  of  5.5"]  since  its  approximate  pole-zero  configuration 

can  be  reconstituted  from  the  tables  once  they  have  been  constructed. 


In  the  next  section  it  is  shown  how  particular  patterns  of  pole-zero  relation¬ 
ships  of  the  booster  relating  function  directly  influence  the  choice  of  the  c^  com¬ 
pensation  unit  for  the  control  system  and  the  closed  loop  performance  of  the  booster 
flight  control  system.  Hence  reference  to  these  tables  aids  in  making  a  quick 
identificati(Mi  of  those  areas  in  which  the  design  of  a  particular  booster  flight 
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Fig.  S(a-d).  Examplts  illustrating  ih*  difftrant  pola-ttro  eanfigurotiont  of  (RF)^ [  ] 

(including  ont  banding  mod#,  and  aswming  ‘ 
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INCREASING  MAGNITUDE  OF  “OPEN  LOOP  GAIN,"  kj,  OF  THE  ROOT  LOCUS  OPERATION 
•.  Peitibl*  lafo  leeotient  for  Com  3  of  Fig.  6(c). 


t - 

- 


—  “*1 


-e- 


-1 

(»)rd 


«- 


g.  A  variation  of  Com  4 


Pig.  5 (•iffS).  Exomplgt  illustrating  tht  dtffsnnt  pols*zsra  configurations  of  (RF)^  [ {.;*.] 
(Including  ono  bonding  mods,  ond  ossuming  \Q|  -  *  ^lg|)* 
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Table  1.  Tabulation  a(  the  phase  (or  range  of  phase)  of  (RF)|,[  ]  /[I  +  p  (SR),^]  evaluated  at  its 

singularities  in  the  upper  holf  of  the  complex  plone.  This  is  for  the  cases  of 
(RF)b[  S;S*.]  illustroted  in  Fig.  5  o-g. 


■ 

Conditkmi 

PhoMjIthe  singulaiitiH 
(indai 

«l<RFli,[j{S*]/[l+p(SR)rt] 
poet  of  la) 

RolHive  VilMof 

Roiaivo  Vilueof 

Ri^Body 

Polei  of  the 
ithBondioiModo 

Hi^  Fiogiacy 
orRHPZdios,M 
ZomsaRed  Axis 

Lowor  Ftoqancy 
or  LHP  Zms 

1 

<  1 

<0 

0® 

90® 

90® 

90® 

2 

<  1 

>  0 

0® 

270* 

90* 

or 

0*  and  IBCP 

270* 

3 

>  1 

<  0 

0® 

270* 

90®,  or 

180®  to  0®,  or 

270* 

270? «» 
0®to-180®,  or 
90® 

■ 

>  1 

>  0 

0® 

_ 

90® 

90* 

or 

0*  and  wf 

90® 

RHP  ■  iHM  had  ImM  of  Ike  CMptex  piM 
LHP  ■  M  hatd  hdl  of  Hie  coalex  piaie 


TobU  2.  Tabulation  of  the  possible  locations  of  the  zeros  of  (RF)|,[  ]  /  d  P  (SR),^] 

in  the  upper  holf  of  the  complex  plone  (or  the  coses  illustrated  in  Fig.  5  o.g. 


CASE 

LOCATION  OF  THE  ZBIOS  ASSOCIATED 
IITH  THE  Hh  MODE 

LOCATIOH  OF  THE  •*TAIL*»AG$-006" 
ZEROS  OF  THE  BOOSTER 

1 

On  the  imaginary  axis  between 
oij  >/Ij  and  0 

On  the  imaginary  axis  between 
%  and  cujj 

2 

On  the  imaginary  axis  between 

ci)g^  and  a>i 

On  the  imaginary  axis  between 

u-  and  °°,or  on  the  real  axis 
*0  ' 

3 

Both  zeroa  are  found  on  the  imaginary  axis  between  euj  and  ai2^;  or 

One  zero  occurs  in  the  RHP  and  the  other  occurs  in  the  LHP  between 
the  greater  of  w;  >/i{  and  q>^^,  and  0;  or 

Both  zeros  are  found  on  the  imaginary  axis  between  and  0. 

■ 

On  the  imaginary  axis  between 
ti>i  Vlj  and  oij. 

On  the  imaginary  axis  between 
and  °^or  on  the  real  axis 

RHP  ■  if|M  had  hrif  of  tiM  coegilix  pim 
UP  •  loR  had  hdf  of  tho  conpla  plae 
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control  system  may  be  troublesome  and  those  areas  where  design  will  be 
straight  forward. 

RELATION  OF  CONTROL  SYSTEM  ANALYSIS  TO  CONTROL  SYSTEM  DESIGN 


The  design  process  consists  of  matching  the  sign  of  X  /0 

g{ 


with  a 


tK 

compensation  which  provides  the  proper  phasing  of  the  i  mode  feedback 

XU 

signal  to  produce  a  maximum  of  artificial  phase  stabilization  for  the  i 
mode.  Alternatively,  the  compensation  may  provide  sufficient  attenuation 
of  the  i^^  mode  feedback  signal  to  insure  stability  of  the  mode  regardless  of 
the  sign  of  the  i^^  bending  mode  pickup  and  the  phase  shifts  due  to  the  remainder 
of  the  system.  This  process  is  carried  out  for  all  the  modes  which  are  included 
in  the  analysis.  A  combination  of  the  above  compensation  methods  may  be 
employed  for  stabilization  of  a  mode  with  success;  however,  such  a  system 
requires  extensive  study  to  insure  the  preservation  of  specified  stability 
margins  of  both  phase  and  gain. 


As  previously  shown  the  locations  of  the  zeros  of  the  open  loop  booster 
relating  function  are  influenced  by  the  relative  values  unity 

*iid  X  .  /  to  zero.  The  effect  of  some  of  these  possible  locations  on 

**gi  gi 

the  system  dynamics  and  general  design  procedure  is  illustrated  by  a  simple 

example  considering  the  rigid  booster  and  one  bending  mode.  The  nozzle  and 

actuator  dynamics  are  represented  as  a  single  first  order  lag;  this  assumes 

that  the  second  order  dynamics  associated  with  the  nozzle-actuator  linkage 

are  negligible.  The  rate  and  rate  integrating  gyros  are  located  at  the  same 

station  along  the  booster  body.  (  X—  =  X  =  X,  )  The  performance  functions 

*^*1  *1 

of  c^  and  c^  compensation  elements  are  equal  to  unity.  The  dynamics  of  the 
forward  path  compensation,  (PF)  ,  are  used  to  phase  compensate  the  system 

used  in  this  example. 


The  performance  function  for  the  closed  loop  control  system  can  be  devel¬ 
oped  from  Figs.  1  and  3.  This  is  given  below. 
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S..  {(Pa(S.;Sl  iff).,} 

(PF)k..tfl.;81  -  ,  *s..  «PF),[,^,8,  (RF)ks;8;i  (PF).,| 


The  characteristic  roots  of  the  booster  flight  control  system  are  the  poles  of 
(PF)  Q  gj.  These  may  be  determined  from  the  above  equation  using 

root  locus  techniques. 

If  the  departure  and  arrival  angles  for  this  root  locus  are  calculated  for 
a  case  without  compensation*  (PF)  =  1,  the  approximate  phase  angle  required 

of  (PF)  at  the  bending  frequency  to  maximize  the  negative  real  coordinate  of 

the  closed  loop  bending  poles  becomes  apparent.  A  suitable  type  of  compensation 
may  then  be  chosen  to  provide  the  needed  phase  shift  at  the  bending  frequency 
without  causing  the  rigid  body  response  to  change  by  any  substantial  amount. 

Figure  6  presents  sketches  of  the  root  locus  for  three  cases  of  the  possible 

locations  of  the  zeros  of  (RF).  ,  ^  with  and  without  second  order  low  pass 

bl  6;  6^ ) 

compensation  to  achieve  phase  shift  at  the  bending  frequency.  It  is  obvious 
that  Figs.  6(a)  ,  (d)  ,  and  (e)  represent  cases  of  correctly  compensated 
systems,  and  that  (b)  ,  (c)  and  (f)  represent  cases  of  incorrectly  compensated 
systems. 

It  is  quite  possible  for  the  booster  flight  control  system  to  change  from 
one  of  the  situations  illustrated  in  Fig.  6  to  another  during  the  flight  of  a  typical 
booster.  For  instance: 

1 .  If  the  bending  mode  resonant  frequency  becomes  greater  than 
the  "tail- wags -dog"  frequency  of  the  mode  during  the  booster 


*  see  Fig.  6(a),  (c),  and  (e) 
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Fig.  <(a.O.  Root  local  diagrams  for  boootor  flight  control  systoms  Including;  boostor  and 
ono  bonding  mods  (sovorol  typoi)  and  control  lystom  with  and  without  compmoation 
(Pp)  bet  [  >1 '  [  »>»;.] 
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flight*(  /  w  ^  <1  changes  to  >1) 

changes  from  the  situation  illustrated  in  Fig.  6(d) 


the  system 
to  that  in  Fig. 


6(f)  . 

2.  If  the  rate  gyro  is  located  such  that  the  mode  slope  sensed  by  the 

gyro  changes  sign  during  the  booster  flight’*'*  (  \  <0  changes 

*‘*1  *1 

to  X  /0  >0)  the  system  changes  from  the  situation  illustrated 

rgi  gj 

in  Fig.  6(a)  to  that  in  Fig.  6(c). 


If  such  a  change  of  booster  parameters  does  occur,  it  is  necessary  to  make 
corresponding  changes  in  the  flight  control  system  parameters  to  preserve 
the  desired  performance. 


A  case  considering  additional  bending  modes  maybe  handled  in  a  manner 
similar  to  the  example.  However,  this  would  impose  additional  constraints 
upon  the  choice  of  (PF)  In  the  case  of  modes  being  phase  stabilized,  the 

approximate  phase  of  (PF)  at  the  bending  mode  resonant  frequencies  would 

*■2 

be  specified  by  the  choices  of  the  signs  for  X^,  and  by  the  reqiairemenk 

that  the  negative  real  coordinate  for  each  of  the  closed  loop  bending  poles  be 
increased.  In  the  case  of  modes  being  amplitude  stabilized,  the  maximum 
allowable  gain  of  (PF)  at  the  bending  mode  resonant  frequencies  would  be 

specified. 


SUMMARY 

This  paper  has  advanced  a  simplified  concept  of  analysis  for  a  flexible 
booster  flight  control  system.  The  main  concern  has  been  the  development  of 
a  simplified,  factored,  open  loop  relating  function  for*  the  booster  control  sys¬ 
tem.  This  provides  the  key  to  the  evaluation  of  system  performance.  It  was 
shown  in  the  closing  section  that  once  this  relating  function  has  been  evaluated, 
the  closed  loop  system  performance  can  be  assessed  using  straightforward  suid 
familiar  techniques. 

*see  first  footnote  on  page  14 

’^’''Changes  in  the  sign  of  a  mode  slope  sensed  by  a  gyro  occur  due  to  the  wide 
variation  of  the  booster  mass  distribution  during  the  flight  profile.  This  is 
very  likely  to  occur  when  an  attempt  is  made  to  place  a  gyro  near  an  antinode 
of  a  particular  mode. 
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This  analysis  Is  unique  in  that  it  enables  this  very  complicated  d)mamics 
problem  to  be  discussed  and  interpreted  Intelligently  in  qualitative  terms.  An 
engineer  thoroughly  familiar  with  this  method  can,  on  a  moments  notice,  be 
prepared  to  discuss  the  control  problems  involved  In  the  preliminary  design 
of  a  particular  booster.  The  knowledge  required  for  discussion  involves  only 
well  known  quantities  such  as  thrust,  moment  of  inertia,  nozzle  mass,  rough 
estimates  of  the  bending  mode  frequencies,  mode  shapes,  and  slopes  for  the 
particular  booster. 

There  are  a  host  of  points  concerning  the  booster  flight  control  system 
problem  which  are  not  suitable  for  presentation  here  with  the  fundamental 
principles,  but  which  may  be  adequately  treated  by  the  methods  of  this  analy¬ 
sis.  Among  these  are  the  effects  of:  elements  of  the  c^  compensation  unit  in 
the  rate  damping  loop  which  do  not  appear  in  the  orientation  control  loop,  ele¬ 
ments  of  the  c^  compensation  unit  in  the  orientation  control  loop  which  do  not 
appear  in  the  rate  damping  loop,  aerodynamics,  propellant  sloshing  modes, 
and  the  couplings  of  bending  and  sloshing  modes  with  the  actuator -nozzle 
motion. 
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APPENDIX  A 

APPROXIMATE  LINEARIZED  EQUATIONS  OP  MOTION 
IN  A  SINGLE  PLANE  FOR  A  FLEXIBLE  BOOSTER 


The  perturbation  equations  of  motion  of  the  booster  in  a  single  plane 
used  as  the  basis  of  this  study  are  essentially  those  developed  in  Reference  2. 
These  have  been  restated  in  a  sli^tly  different  coordinate  system  and  nota¬ 
tion  (standard  aircraft  body  axis  coordinate  system)  in  Equation  Summary  A-1, 
The  coordinate  system  definition  and  the  definition  of  important  system  phys¬ 
ical  constants  are  shown  in  Fig.  A-1. 

These  equations  can  be  further  simplified  while  still  retaining  the  essen¬ 
tial  dynamic  features  resulting  from  the  structural  degrees  of  freedom  of  the 
physical  system. 

The  simplifying  assumptions  are  as  follows: 

1.  Aerodynamics  sure  neglected. 

2.  Cross  couplings  among  the  bending  modes  and  couplings  of 
the  bending  modes  with  the  rigid  body  modes  due  to  engine 
thrust  are  negligible. 

3.  Couplings  of  the  bending  modes  with  the  nozzle  dynamics 
su'e  neglected. 

The  advantages  of  using  these  assumptions  for  simplifying  the  equations 
of  motion  are  as  follows: 

1.  The  rigid  body  mode  and  the  individual  bending  modes  of 
the  booster  are  decoupled  from  one  another  in  the  approx¬ 
imate  equations  of  motion.  This  isolates  the  important 
bending  effects . 

2.  The  relations  between  6  and  the  4  <ind  q^  variables 

(or  their  time  derivatives)  become  simple  second  order  dif¬ 
ferential  equations.  The  approximate  response  of  any  of 
these  variables  may  be  examined  by  cascading  a  simple 
relating  function  with  that  describing  6  for  the  closed  loop 
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booster  flight  control  system.  The  low  frequency  response 
of  these  variables  is  well  approximated  only  when  the  dynamic 
pressure  is  low.  The  high  frequency  response  (the  major 
concern  of  this  report)  is  unaffected  by  neglecting  the  aero¬ 
dynamics. 

3.  When  the  nozzle  d}mamics  are  considered  separately  from  the 
booster  body  motion,  the  effects  of  inertial  reaction  torques 
which  couple  the  rigid  body  and  bending  modes  with  the  nozzle 
dynamics  may  be  described  by  a  feedback  of  each  of  these 
simple,  second  order,  uncoupled  modes  around  the  second 
order  nozzle  dynamics. 

The  introduction  of  these  assumptions  after  elimination  of  the  intermediate 
variables  v,  w,  u^,  and  by  substitution  reduces  the  equations  of  motion 
to  those  of  Equation  Summary  A-2. 


EQUATION  SUMMARY  A-1 
FLEXIBLE  BOOSTER  EQUATIONS  OF  MOTION 


Normal  force  eqaetion 

mo  -  [C.  qS]  a  -  [m  g  cos  0  4  [(C„  +  Cp)  qS]w  -  -  [TlA  (Al-1) 

Components  of  normal  acceleration 

-  [Vp  +  0*1  ®  -  [Vp]  9  (AI-2) 

External  moment  equation 

[Ip2]fl  - -[C„  (p  qS]a4  [Co  qS]v -[C„  Ip  qS]w  +  [710^41 -  [Ti,]« 

“  “  (Al-3) 

Deflection  of  center  line  due  to  nozzle  deflection 

Imlv  - -[|„  m„]a  (Al-4) 

Rotation  of  center  line  doe  to  nozzle  deflection 

(llw  -  -[Igla  (Al-5) 

Deflection  at  the  engine  gimbal  point 

Uj  -  [l]v  +  [|p]w  +  (0g,]q)  +  [^g2l<)2  ■*■•••  +  l^g|lfl|  (Al-6) 

Slope  at  Che  engine  gimbal  point 

0g  -  [l]w  -  [Ap^lq,  -  Ug^lqg  -  - - (Al-7) 

i*!*  bending  mode  generalized  coordinate 

[m  2m  ^1  •*1  p  4  m  «p]  q|  •  -[(T  4Cd  qS)X,|]v-rr  ^gjl^^g-Cinn  <1  P’  +  T^g,]« 

(AI-8) 

Nozzle  moment  equation 

n„  P*  4  2  I„  ^  p  4 1„  «2ia  -  IJ«,  -  ll„  m„]0,  -  ng  P»  4  e„  m„  g  eoa  p]9 
4  IT  4  Cp  qS]  V  -  [|„  m„  p*lMp  ~  D„  P*  +  (m„  l„/m)  (T  4  Cp  qS)]0,  (Al-9) 
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« 


Motion  sensed  by  the  (  )  gyro 

d(  )  -  [p]  +  tp]  w  -  S  [p  A(  )j  1  q| 
i 

-  Ip]5-  1  [pA(  )|]  q| 


(AI-10) 


NOTES;  These  stations  are  modified  from  those  of  Ref.  2  by  a  change  in  the  definition  of  the 
coordinate  system  and  a  change  in  the  number  of  nosdes  swiveled  in  a  single  plane. 
The  effect  of  a  distributed  aerodynamic  normal  force  on  the  bending  has  been  neglected. 
Fuel  sloshing  effects  have  been  neglected. 

In  general,  C„  and  Cq  are  negative. 
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C.O.  OF  TOTAL  BOOSTER 


EQUATION  SUMMARY  A-2 

SIMPLIFIED  FLEXIBLE  BOOSTER  EQUATIONS  OF  MOTION 
Nomal  force  equation 

[m]a,  .  -[mgcos/910-  [t  ^1-  p)]*  (A2-1) 


Components  of  normal  acceleration 

o,  -  (Vp  +  o^la-  [Vpltf 


(A2-2) 


External  moment  equation 


np>i.  .-[n,(u-^)]. 

Generalized  coordinate  of  the  i'^  bending  mode 

(p*  +  2  Cj  <uj  p  +  w*)  -  T Agj  j  q|  - 

Nozzle  moment  equation 

Motion  sensed  by  the  (  )  gyro 

9()  -  <i-[py]«-  S[p^(  q, 

-  e  ~  ),^  'll 


(A2-3) 


(A2-4) 


(A2-5) 


(A2-6) 
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APPENDIX  B 

OBTAINING  THE  FLEXIBLE  BOOSTER  RELATING 
FUNCTION  WHICH  INCLUDES  COMPENSATION 


9* 


The  method  used  in  Equation  Summary  1  to  develop  the  relating  function 
(RF)|^^  6  -6"]  been  extended  to  include  different  compensation  elements  in 

each  of  the  gyro  signal  paths.  The  results  are  summarized  in  Equaticm  Sum¬ 
mary  B-1. 


The  compensation  used  in  each  gyro  path  is  required  to  be  in  factored 
form  and  the  numerator  and  denominator  orders  must  be  equal.  If  the  number 
of  finite  poles  should  exceed  the  number  of  finite  zeros  in  the  actual  compen¬ 
sation  being  investigated,  additional  zeros  must  be  carried  at  infinity  to  equal¬ 
ize  the  numerator  and  denominator  orders.  The  rate  path  compensation  is 
given  as  (PF)^  »Nj(p)/Dj(p),  where  the  order  of  N^  and  D^  is  p'*;  and  the  rate 

integrating  path  compensation  is  given  as  (PF)^  >  Ng(p)/Dj(p),  where  the 

V  ^  . 

order  of  Ng  and  Dg  is  p  .  The  number  of  zeros  associated  with  6-6"] 

when  written  as  a  raticmal  polynomial  is  (2n  +  u  v  +  1  2).  Some  of  these 

zeros  may  occur  at  infinity  as  previously  mentioned.  The  (2n  u  +  v)  zeros 
result  from  placing  the  polynomials  over  a  common  denomixiator;  two  "tail- 
wags-dog"  zeros  from  the  inertial  reactions  of  the  gimballed  nozzle;  and  one 
zero  from  the  rate  gyro  "differentiation'.  In  addition  it  is  required  that 
Nj(0)/Dj(0)  »  Ng(0)  /Dg(0)  -  1. 


This  additional  generalization  increases  the  complexity  of  the  expression 
for  (RF)^I  6-6"]  substantially  because  of  factorizations  of  (^3^2  ^  [const.] 

pN^Ug)  which  must  be  performed  for  each  mode  included  in  the  analysis.  (See 
Eqs.  Bl-1,  Bl-2  and  the  associated  auxiliary  equations  of  Equation  Summary 
B-1.) 
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The  equations  (RF)b(j.  with  (PF)^^  ^  1  presented  below 

without  derivation.  The  method  of  derivation  is  essentially  similar  to  that  of  Equation 
Summary  1. 

The  relating  function  including  the  1st  bending  mode  in  parallel  with  the  rigid  body  is 


(R^blS:  S'] 


H[S;S1 _ 

i.  2^,  p2 


\  O^tui  Oitjf/I  1 


B,  1 

Xi  -r 

’  A,  1 


T?r5j’ 


(Bl-1) 


A,{p)  -  (N3  D,  >  p  (SR),j  N,  D3)  fw  p  +  ^  (Bl-la) 

\  0^  ftj,  0^  ft)'  /  \  / 


B,(p)  .  p»^N,D,.p(»)rf^N,D,) 


(BMb) 


“^Ipi  ^b[S;qi] 


‘b[S;6il 


(Bl-lc) 


Equation  Summary  B*1  (Page  1  of  2) 

Equations  for  the  Flexible  Booster  Relating 
Function  Which  Includes  Compensation 


The  reletiag  function  for  n  bending  modes  can  be  obtained  by  performing  the  following 
operations  on  the  relating  function  which  includes  (n-1)  bending  modes  (n  >  2,  3,  . . . ). 
The  effects  of  the  bending  mode  will  be  added. 


(RF)b  [«;«']  • 


^hlB;  9] 


'  \  Oi  Oi  fti?/ 


p2D,  Dj  n(l  p+; 

'  Oj  Oj 


Kjji  '' 


-I 


(Bl-2) 


where 


.  .  (  the  (2n  -f  u  4-  y  4  1)  zeros  of  Wl 

'  **  “  for  (n-1)  bending  modes/  ^  ^  “  j  *>j  ®  j 


B,(p,  .  pP  (n.  D,  .  p  ,SR,„  ^  N.  D,)  (,  .^)"n’  ('  -^P  * 


(Bl-2b) 


Kj  - 


-A| 


«i 


T 


b(d:  9] 


(B1-2C) 


The2n  -t  u  -t  v  1  zeros  of  Eq.(Bl*2s)must  be  expressed  in  the  form 
(1  +  rp). 


( 


l+lip  + 
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Equation  Summary  B-1  (Page  2  of  2) 

Equations  for  the  Flexible  Booster  Relating 
Function  Which  Includes  Compensation 


GLOSSARY 

SYMBOLS 


^0 

D 

I 

Ie 

Jn 

Ki 

N 

q„,] 

(RP)  r-  .  1 

"W|ni 

s 

;  a,] 
(SR)rd 


T 

V 

"I 


•etodynatnic  drag  coefficient,  Co 


aerodynamic  normal  force  coefficient, 
aerodynamic  drag  force 


dN/da 

qS 


total  moment  of  inertia  of  booster  about  the  center  of  mass 
Jfi  +  "n  *9 

moment  of  inertia  of  nozzle  about  gimbal  point 

non-dimensional  **open  loop  gain’*  used  in  the  root  locus  factorization  of  the  numerator 
of  (RF)^  [  j.j*  ]  when  the  effects  of  the  ith  bending  mode  are  being  added. 

aerodynamic  normal  force 

performance  function  of  component  n,  relating  the  measurable  output,  to  the 
measurable  input,  q|„ 

relating  function  describing  the  mathematical  relationship  between  the  output, 
and  input,  q|„,  of  component  n 

reference  area 

sensitivity  of  the  orientation  control  loop  (5,^  ) 


sensitivity  of  the  rate  damping  loop  (S^j) 
rate  damping  sensitivity  ratio  (SR),^  > 


static  sensitivity  of  the  performance  function  or  relating  function  of  component  n 

relating  the  output,  to  the  input,  q|„ 

thrust 

booster  velocity 

modification  factor  for  natural  bending  characteristics  of  the  mode  when 
the  rocket  engine  is  thrusting  Oj  -  1  -  '^^*1 
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“x,  y,  or  I 
9 

•‘I 

^0 

<n 


m 

•"i. 

P 

q 

<11 

t 

"( ) 

V 


p 

8 


«a 

«e 

^  ) 
$ 

e 


■  acctleratioo  of  booster  center  of  mass  in  the  positive  x,  y,  or  i  direction 
••  acceleration  of  gravity 

-  dimensional  form  of  K|;  the  form  of  the  open  loop  gain  actually  used  in  performing  root 
locus  factorizations. 
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distance  between  the  booster  center  of  mass  and  the  nozzle  gimbal  point 

distance  between  the  center  of  mass  of  the  nozzle  and  the  gimbal  point 

distance  between  center  of  pressure  and  the  center  of  mass  (positive  for  center  of  ptesaure 
forward  of  the  center  of  mass) 

total  mass  of  booster 

mass  of  the  nozzle 

Laplace  operator;  a  complex  number  p  ■  o  iu  (o  and  a>  ate  real  numbers,  j  *^1  ) 
dynamic  pressure 

generalized  deflection  coordinate  of  the  bending  mode 
time 

deflection  from  the  uodefoimed  elastic  axis  at  the  station  identified  by  the  subscript 

translation  of  the  rigid  body  from  the  i  axis  in  the  z  direction  due  to  negative  nozzle 
deflection  (T  ■  0) 

rotation  of  the  rigid  body  centerline  from  x  axis  due  to  negative  nozzle  deflection  (T  -  0) 
angle  of  attack  -  angle  from  the  velocity  vector  to  the  undeformed  elastic  axis 
angle  from  the  reference  trajectory  to  the  local  vertical 

nozzle  deflection  angle  from  the  centerline  of  booster  at  the  gimbal  point  to  the  centerline 
of  the  rocket  nozzle 

hydraulic  actuator  output  in  units  of  equivalent  nozzle  deflection 
commanded  nozzle  deflection 

damping  ratio  or  effective  damping  ratio  of  the  mode  indicated  by  the  subscript 

angle  from  the  reference  trajectory  to  the  undefotmed  elastic  axis 

angle  from  the  reference  trajectory  to  the  rigid  body  centerline  (nozzle  deflected) 
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fig  >  coomaiuled  ffngle  from  refetMce  trajectory  to  undeflected  booster  centerline 

a 

0(  )  ■>  an|ular  velocity  sensed  by  s  gyro  located  at  station  indicated  by  subscript 

-  normalized  bending  mode  slope  amplitude  at  station  identified  by  subscript 

).-  normalized  i'^  bending  mode  deflection  amplitude  at  station  identified  by  subscript 

d'(  )  ■  slope  of  the  booster  centerline  with  respect  to  the  undeformed  elastic  szis  at  the  station  in¬ 
dicated  by  the  snhscript 

ai(  )  >  undamped  naturai  frequency  or  effective  undamped  natural  frequency  of  the  mode  indicated  by 
the  subscript 


S' 


(  ) 


NOTES:  (1)  Station  numbers  are  distances  from  a  reference  point  on  the  undeflected  center  line.  This 
rcfetence  is  usually  chosen  forward  of  the  nose  and  is  positive  in  the  opposite  sense  to 
the  X  axis. 

(2)  All  angles  are  defined  as  positive  for  a  small  angle  vector  treatment  in  which  the  vector 
sense  is  positive. 


SUBSCRIPTS 

0  ■  control  system  actuator  or  servo 
i  •>  i***  bending  mode  (I  ->  1,  2,  3,  •  •  •) 
ig  >  rate-integrating  gyro 
b  m  booster 

bcs  >  booster  flight  control  system  inclnding  booster 
n  a  nozzle 
rg  a  rate  gyro 

*0t  Z|  a  “tail-wags-dog”  zero  of  rigid  body  or  i^^  bending  mode 
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